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A new approach to the proof of the Arrhenius formula of kinetic theory is proposed. We 
prove this formula starting from the equation of diffusion in a potential. We put this diffusion 
equation in the form of evolutionary equation generated by some Schrodinger operator. We 
show that the Arrhenius formula for the rate of over the barrier transitions follows from the 
formula for the rate of quantum tunnel transitions for the considered Schrodinger operator. 

Relation of the proposed approach and the Witten method of the proof of the Morse in- 
equalities is discussed. In our approach the Witten spectral asymptotics takes the form of the 
low temperature limit and the Arrhenius formula is a correction to the Witten asymptotics. 



> 
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1 Introduction 

The Arrhenius formula [T] for the reaction rate in chemical kinetics describes chemical reactions 
which are related to the over barrier transitions between the potential wells. With this formula 
the reaction rate is proportional to 



where AE is the difference of the energies of the transition state (the activation barrier) and of 
the initial state. 

The Eyring formula for the reaction rate (the modification of the Arrhenius formula) 

replaces the difference of energies by the difference of the free energies AE. 
Let us remind that the free energy F{G) of the set of states G is defined as: 

e -0F(G) = f e -m*) dx . 

Jg 

Here U(x) is the energy of the state x. 
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In the present paper we discuss the proof of the Arrhenius and the Eyring formulas starting 
from the equation of diffusion in the potential 



-± = Af + pVf-VU + PfAU, 



(2) 



where x G M d , / = f(x,t) is the distribution function, U = U(x) is the potential (a real valued 
function), U G C 2 (M. d ), (3 = 1/kT > is the inverse temperature. 

We show that the above equation is equivalent to the evolution equation generated by some 
Schrodinger operator. The ground state for this operator corresponds to the Gibbs state for 
the diffusion in the potential (|2J). The energy of the first excited state describes the velocity of 
relaxation to the Gibbs state. 

For the case of the double well potential the energy of the first excited state can be estimated 
quasiclassically. In the quasiclassical approach the splitting of the ground state for the Schrodinger 
operator is generated by the tunnel transitions between the potential wells. We show that this 
estimate for the energy of the first excited state gives the Arrhenius formula. 

Therefore we prove the Arrhenius formula which describes the over barrier transitions in the 
classical kinetic theory using the spectral properties of some Schrodinger operator related to 
quantum tunneling under barrier transitions. 

The relation of the Arrhenius formula and non-Arrhenius relaxation on complex landscapes 
of energy was discussed in many works, in particular, see pQ. In [2] the relation of dynamics on 
complex landscapes and p-adic diffusion was considered. In [3J, jl] the quantum low temperature 
limit for the reaction rates which is related to tunneling transitions was discussed, see also [5] 
for the discussion of the quantum corrections to the Arrhenius formula. In the present paper we 
use quantum tunneling for the corresponding Schrodinger operator to prove the purely classical 
Arrhenius formula. 

The structure of the present paper is as follows. 

In section 2 we discuss the diffusion equation in the potential and construct the corresponding 
Schrodinger operator. We also prove the stabilization property for the evolution described by the 
considered diffusion equation. 

In section 3 the quasiclassical approach to the introduced in section 2 Schrodinger operator 
is discussed. We show that in the regime of low temperatures the bound states correspond to 
potential wells related to local minima of the potential U (i.e. the potential in the diffusion 
equation (J2J)). We prove that the quasiclassical formula for the splitting of the ground state for 
the double well potential in the low temperature limit implies the Arrhenius and Eyring formulas. 

In section 4 we prove the Eyring formula beyond the frameworks of the quasiclassical approx- 
imation. 

In section 5 we discuss the relation of the approach of the present paper and the Witten 
approach in the Morse theory. 

2 The diffusion equation and the Schrodinger operator 

Equation (j2J) can be put in the form 



dt 



div [e-^grad [fe pu ]]- 
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The equilibrium state fix) = const e _/3C/<x ) (the Gibbs distribution) is the stationary solution of 
the above equation. 

With the substitution ip = e^ u ^ 2 f this equation is put in the form 

= -Hi(>, H=- [eW 2 Ve-^ /2 ] • [e-^ 2 Ve^ 2 ] . (3) 
The operator H is self-adjoint and positive, since H is a finite sum of positive operators 

H = A* ■ A = A*A U A = e-^VeW A . = e -wd^. (4) 
i=i * 

The gradient part of H cancels and we get the Schrodinger operator 

H = -A + V, V = -^AU + ^{VU) 2 . (5) 

Here V is the effective potential. The spectrum of the Schrodinger operator H is bounded from 
below by zero. Since the function e~^ u l 2 satisfies 

He-^l 2 = 0, 

then, if e~^ u l 2 e L 2 (M d ), i.e. the Gibbs state is integrable e~^ u E L 1 (R d ), zero will be the minimal 
point of the discrete spectrum of H and e~^ u ^ 2 will be the ground state. 

Let also the operator H satisfies the condition of the existence of the energy gap: the spectrum 
of H is a subset of {0} (Jt-^i) 00), £1 > (we assume that E\ belongs to the spectrum). The value 
E\ we call the energy gap for the operator H. 

The positivity of H and the existence of the energy gap implies the convergence of the solution 
of equation ([2]) to the Gibbs state with t — > 00. We formulate the following theorem (a variant of 
a theorem about stabilization of solutions of parabolic equations). 

Theorem 1 Let for equation (Tj2)) the potential satisfies U G C 2 (M. d ) and the Gibbs state is inte- 
grable 

r e-^^dx < 00. 



Let also the potential U satisfies the condition of the existence of energy gap. 

Consider for equation (0) a Cauchy problem in the space L 2 (R d ,e^ u ^dx) of quadratically 
integrable with the weight e" u functions. The initial condition fo for a Cauchy problem belongs 
to this space, i.e. fo satisfies 

r \f (x)\ 2 e mx) dx < 00. 



Then the solution f(t,x) of the described Cauchy problem for equation (TJJ) with the initial 
condition fo(x) exists, is unique and with t — > 00 tends to the function proportional to the Gibbs 
state 

ton f(x,t) = e-M*) V'Wte 



t— >oo 
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Proof The scalar product of the functions ip(x,t) = f(x,t)e@ u ( x '/ 2 and ipo{%) = e P u ( x >' 2 in 
L 2 (R d ) 



is 



I? 



f(x, t)dx. 



This implies that the solution / of ([2]) satisfies the condition of conservation of the number of 
particles 

(/ f f(x,t)dx = (6) 



dt 

because the energy of the ground state ip is equal to zero. 

Then since the energy of the ground state is zero and there exists the energy gap 



t— >oo 



t— ¥00 



lim f(x,t) = lim i;(x,t)e-^ 2 = \im (i;^ ) LHRd) M^e~ mx)/2 
By §6§ this implies 



lim f(x, t) = e 



t— >oo 



-Mx) I fo( x ) dx 



which finishes the proof of the theorem. □ 



Let us discuss the behavior of the solution for the case when the operator H has the discrete 
spectrum and the initial condition ip(0) for ([3]) belongs to L 2 (M. d ). Let the initial condition 
possesses the expansion ^(0) = ^ over ^ ne eigenvectors of H, which correspond to the 

eigenvalues Ei < E i+ i, ^ is proportional to e~^ u ^ 2 , E — 0. Then 



i=2 



Therefore the evolution %fj{t) tends to the ground state ipo proportional to e _/3C///2 , and the rate of 
the convergence for large t is described by the term ^ 1 e~ £ ' lt in the above expansion corresponding 
to the first excited level of the Hamiltonian H. 



Example Let us discuss the case of the quadratic potential U = a\x\ 
Schrodinger potential V is also quadratic: 



Then the effective 



V 



^AU + ^(VU) 2 



-(3ad + (3 2 a 2 \x\ 2 . 



Here d is the dimension of the space. 



d 

H = -A + V = J2' 2 ^- 



(3a dx\ 



+ (3 'ax 2 



We have d harmonic oscillators with the eigenvalues Ei = 2(3arii, rii = 0, 1, 2, In particular 

the term — §AC7 cancels the energy of the vacuum and makes the energy of the ground state for 
H equal to zero. 

The energy gap (between the ground and the first excited states) in this model is equal to 2f3a. 
In the considered single well case the Arrhenius formula is not applicable. 
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Remark 2 We consider the effective potential V = — f A£7 + ^-(V£/) 2 in the regime of low 
temperatures (large f3) and a generic potential U . Generic here means that U is a Morse function 
i.e. it possesses a finite number of critical points (the points where Vf/ = 0) and the Hessian (the 
matrix of second derivatives of U) in the critical points is non-degenerate. 

Since in the low temperature regime the leading contribution to the effective potential comes 
from the term containing the square of the gradient of U, the minima of the effective potential V 
will correspond to the critical points of the potential U. In this case the values of the effective 
potential V in the minima corresponding to the minima of U will be negative (by the positivity 
of the second derivatives in the minima of U). In the minima of V corresponding to the maxima 
of U the values of V will be positive. For saddle points of U the values of V in the corresponding 
minima will depend on the sign of the Laplacian of U. 

In the next section we show that the minima of the potential U satisfy the Bohr-Sommerfeld 
condition of existence of bound states for V . 



3 Quasiclassical approach 

In the present section we consider the quasiclassical approximation for the Schrodinger operator 
under consideration (in this section we restrict ourselves to the one dimensional case and the 
double well potential). We show that in the low temperature regime in this approximation there 
exist states of the discrete spectrum for the wells of the effective potential V which correspond 
to the local minima of U and the splitting of energy for the ground state related to transitions 
between these potential wells is described by the Arrhenius formula. 

Remark 3 The quasiclassical Bohr-Sommerfeld quantization condition has the form [6] 

If 1 

p(x)dx — n + -, 



2nh J y ' T 

where the integration runs over the periodic orbit of the quasiclassical particle, n — 0, 1, 2, . . . is 
the number of the quantum state. 

Let us investigate, using this condition of possibility of existence of a bound quantum state in 
the potential well for operator (JSJ). For the ground state we have the energy E = 0, the momentum 
of the quantum particle \p\ = y/2m(E — V) = \f—2mV with the normalization h = 2m = 1. 

We assume that the bound state in the potential well exists if the Bohr-Sommerfeld condition 
can be satisfied for n — 0, i.e. 

- [ ^V{x)dx > i (7) 

TT J-V>0 2 

where the integration runs over the domain where the potential V = —^AU + ^-(VU) 2 is non- 
positive. 

Minima of the potential V are described in remark [2j For large (3 (low temperatures) minima 
of V correspond to critical points of U. Minima of V related to maxima of U do not satisfy the 
condition since for these minima the potential V is positive. 

For minima of V corresponding to the minima of U the potential V is negative. We get in the 
one-dimensional case, expanding the potential U in the vicinity of the minimum x 



' / !(3 U"(x )-^-(U"(x )x) 2 dx = -[ Vl^dx = \. (8) 



7T ./_V>o V 2 4 7T J_ x 2 
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Therefore for these minima of V the condition ((Tj) will be satisfied. 

In the following we will ignore the potential wells of V corresponding to the critical points of 
U which are not local minima. 



Let us consider the one dimensional double well potential which contains two symmetric po- 
tential wells separated by the potential barrier. Tunnel transitions between the potential wells 
imply the splitting of the energy of the ground state i.e. the ground state becomes the pair of 
close energy levels and the difference of the corresponding energies is related to the rate of tunnel 
transitions between the potential wells. 

For the dynamics generated by the Schrodinger operator with this potential the rate of relax- 
ation to the ground state for large times will be described by the energy Ex of the first excited 
state i.e. by the splitting of energy of the ground state due to the tunnel transitions. Let us show 
that for the Schrodinger operator under consideration this energy splitting will be described by 
the Arrhenius formula. 

The formula for the splitting of the energy of the ground state in the symmetric one dimensional 
double well potential in the quasiclassical approximation [6], Chapter 7, section 50, problem 3, 
has the form 



Ex - E = exp (~ ^ \p{x)\dx\ 



Here p(0) is the momentum of the particle at the barrier of energy and Vol is the volume of 
one of the potential wells. The integration runs over the interval of the tunnel transition where 
the potential V is positive i.e. | U' 2 > U" . The momentum of the quantum particle \p\ = 
^2m{V — Eq) = y/2mV with the normalization used h = 2m = 1. Since Eq = we have 
Ei — E — Ex. 
We get 

El jm exp (_£ ^) jm exp (_/^ 

We are interested in the low temperature limit when the transitions between the potential 
wells are slow. In this case the inverse temperature /3 is large and we can omit the first term under 
the square root in the expression for E\ and obtain the expression 

exp [ — ^ I \U'(x)\dx 
\ 2 Jiu ,2 >u" , 



Vol \ s je.v>*>u" 
Integration runs over the symmetric domain of the potential barrier, therefore 



\U'(x)\dx = 2 / U'(x)dx = 2AU, 

U' 2 >U" J §U' 2 >U";x<0 



where the energy barrier AU is equal to the difference of values of U at the potential barrier (for 
x = 0) and in the well (at the point where ^U' 2 = U"). Since the inverse temperature /3 is large 
the condition ~ U' 2 = U" is approximately satisfied at the critical point U' — 0. 
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This implies that 



Bl = M exp( _ MC/) 



i.e. the energy splitting of the ground state for the operator (jSJ) and correspondingly the relaxation 
rate for equation (J2J) possess the Arrhenius dependence. 

Taking into account the multiplier |p(0)|/Vol we can get the Eyring formula. We have 



l/'(0)| = \l~U"(x ) 



where Xq is at the energy barrier. Modeling the potential in the vicinity of the energy barrier by 
the Gaussian expression AU exp(— a 2 (x — x ) 2 ) we get |p| = a-^fiAU and 

Ei = v / ^AZ7^exp(-/3Af/). 

Since In a is the entropy of the transition state and In Vol is the entropy of the initial state (the 
potential well), the above formula for the energy of the first excited state for the double well 
potential coincides with the Eyring formula (one can ignore the slow dependence on y//3AU). 



4 The Eyring formula 

In the present section we study the splitting of the energy levels for the Schrodinger operator (JSJ) 
in the case of the double well potential U in the multidimensional case. Since we know the exact 
expression for the ground state we can perform the computations analogous to the computations 
made in the previous section beyond the frameworks of the quasiclassical approximation. 
Let us consider the normalized ground state (with the norm one) 



4u(x) 



[J e-M*)dx] 



and the normalized first excited state ipi for the Schrodinger operator (jSJ). These states satisfy 
the equations 



HI* = (-A - ^AU + ^(Vf/) 2 ) = 0, (9) 

Hik = (-A - ^AU + ^(W) 2 ^) Vi = E^. (10) 

Let us multiply the first equation by ipi and the second equation by ipo and subtract the second 
equation from the first one: 

i)iAip Q - V'oA^i = EiipxipQ. 

Let us construct a surface S which goes through the saddle point (the transition state be- 
tween the two potential wells) and separates all the space in the two domains which contain the 
corresponding potential wells of U. Namely the surface S goes through the saddle point and the 
gradient VU is tangent to S (equivalently is orthogonal to the normal to S)0. 



^or the dimension of the space larger than two such a surface is not unique. We choose some surface of this 
kind. 
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We integrate over one of the domains G separated by S (say over the domain which contains 
the first potential well). Applying the Green formula we get 

(-0iA-0o - ipoAipi)dx = I (^ 1 ^^-^ ^-jdS = E 1 j ipxip dx. 

G j QG \ / */ G 



Since ipo is proportional to e 2 u the derivative of -ipo with respect to the normal to S is equal to 
zero: ^ = 0. 

an 

We get the following (exact) expression of the energy of the first excited state using the wave 
functions of the ground and the first excited states 

J G ^iijj dx 

Expression ( FlTT) for the energy of the first excited state contains the integral in the denominator 
which can be estimated as 

Modx= [ $dx= [ e- pu{x) dx = e~ m \ 
g Jg Jg 

since the ground and the first excited states in one of the potential wells can be considered 
approximately equal. Here Fq is the free energy of the domain G (the main contribution to the 
free energy for large /3 comes from the potential well). 

Let us discuss the form of the state We make in ( flOl) the substitution ip\(x) = 9(x)ipo(x). 
The function 9 should be close to constants in the potential wells and should have a large gradient 
in the transition states — barriers between the potential wells. We get (taking into account that 
ipo is proportional to e~ l3U ^ 2 ) the equation 

-AO + (3V9 ■ W = Ei9. 

In the one dimensional case this equation reduces to 

-9" + 09'U' = Et9, 

where E\ is a small parameter. Omitting the small right hand side of the above expression we get 
the approximate solution 

e'( x ) = Ce mx) 

We choose the normalization constant C using the following statement: the function 9 changes 
from one to minus one when moving from one potential well to the other. This implies 

2e PU(x) 



9'(x) 



where the integration runs from one potential well to the other. In the following we will use the 
above approximation for the multidimensional case. 

The integral over the surface S = dG in the numerator can be estimated as follows. The 
expression under the integral is approximated by 

. , ^i(x) of . . ;2 . . 2ef> u M , 2 *> 



On v-;rov-y j e pu(x) dx ^o f e fiu(x) dx 
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The sign changes since the normal to S is directed to the opposite direction. In the following we 
ignore the coefficient two in the above expression (since this coefficient is beyond our approxima- 
tion) . 

For the transition state the integral in the numerator of (llip is supported in the vicinity of the 
energy barrier (where the function 6' is supported). We get the approximation 



-[ -l 



dS ~ / / e- f>u WdxdS = e- fiFL , 

dG 



where F\ is the free energy of the transition state. 
The above discussion implies the following: 

Proposition 4 The splitting E% of the ground state for the Schrodinger operator (03) for the 
double well potential is approximated by the Eyring formula 

E 1 = exp(-P(F 1 -F )). 

Here Fo is the free energy of the potential well (the initial state of the corresponding kinetics) and 
F\ is the free energy of the transition state. 

We proved (at the physical level of rigour) that the transition rate for the kinetics related to 
the diffusion equation (j2J) is described by the Eyring formula. In our approach the rate of over 
barrier transitions is described by the spectral theory of the corresponding Schrodinger operator 
which itself is related to the rate of quantum tunneling transitions between the potential wells. 



5 The Witten Laplacian and the Morse theory 

The Witten approach to the proof of the Morse inequalities [7], see also [8], [9], Chapter 11, is 
related to computation of the spectral asymptotics of the deformed Witten Laplacian L t , 

Lt = d* t d t + dtd*, dt = e~ tf de tf , d* t = e tf d*e~ tf . 

Here t 6 M, / is some Morse function, i.e. a smooth function with the finite number of critical 
points (points where the gradient is equal to zero) where the Hessian (the form of second deriva- 
tives) is non-degenerate, d is the exterior derivative (acting on differential forms), d* is the Hodge 
conjugate to d, t £ K. 

The asymptotics of the spectrum of the deformed Laplacian is computed quasiclassically in 
the limit t — > oo. The corresponding estimates for the Betti numbers coincide with the Morse 
inequalities. 

The Witten approach is connected to the discussed in the present paper spectral properties of 
the Schrodinger operator in the following way. The Witten asymptotics is the low temperature 
asymptotics and the Arrhenius formula is a correction to the Witten spectral asymptotics. 

The operator A of the form fl4]) is a restriction of the deformed exterior derivative dt to functions 
(i.e. zero-forms). We have / = U, t = (3/2. The operator (jSJ takes the form H = d*dt of one 
of the contributions to the deformed Laplacian, the asymptotics t — > oo is the low temperature 
limit. 
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The Witten approach discusses the spectra of the oscillators related to all critical point of /. 
In the present paper we consider only the minima of U (not all critical points). The discussed 
here splitting of the ground state due to tunnel transitions between the potential wells (described 
by the Arrhenius formula) is a correction to the Witten spectral asymptotics (and vanishes in the 
proof of the Morse inequalities). 
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